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Effects of Cosmic String Velocities and the Origin of Globular Clusters 


Ling Linj^ Shoma Yamanouchij^ and Robert Brandenbergei[^ 

Department of Physics, McGill University, Montreal, QC, H3A 2TS, Canada 

With the hypothesis that cosmic string loops act as seeds for globular clusters in mind, we study 
the role that velocities of these strings will play in determining the mass distribution of globular 
clusters. Loops with high enough velocities will not form compact and roughly spherical objects 
and can hence not be the seeds for globular clusters. We compute the expected number density and 
mass function of globular clusters as a function of both the string tension and the peak loop velocity, 
and compare the results with the observational data on the mass distribution of globular clusters 
in our Milky Way. We determine the critical peak string loop velocity above which the agreement 
between the string loop model for the origin of globular clusters (neglecting loop velocities) and 
observational data is lost. 


I. INTRODUCTION 

We have recently [T] made the hypothesis that string 
loops arising from a scaling network of cosmic strings 
seed the formation of the globular clusters which are 
observed to be distributed in the halos of galaxies, in 
particular our own Milky Way galaxy. Our model easily 
explains the observational facts that globular clusters 
are the oldest and most compact star clusters in our 
galaxy, and that they are distributed throughout the 
halo as opposed to only in the disk. 

Cosmic string loops arise via the interaction of infi¬ 
nite string segments which in turn are generated during 
a symmetry breaking phase transition in the early uni¬ 
verse. Cosmic strings are predicted in a large class of 
models of particle physics beyond the Standard Model. 
According to the cosmic string scaling solution [2], the 
distribution of infinite string segments is independent 
of the cosmic string tension fi (which is, in the natu¬ 
ral units we use, equal to the mass per unit length). 
As a consequence, cosmic string loops are also formed 
with a number density which is independent of pL. The 
distribution of string loops is determined by the string 
scaling solution, and depends on /r only through the de¬ 
pendence on /X of a critical loop radius Rc below which 
the number distribution of strings becomes constant [5] . 

In our previous paper ^ we have shown that if we fix 
the one free parameter in our theoretical model, namely 
the string tension yi, then the peak number density and 
the mass distribution are fixed. Demanding that the 
mass distribution peaks at a value corresponding to the 
peak in the observed mass function of globular clusters 
in the Milky Way gave us a value of G/i ~ 10“®'^ (where 
we - as is standard in the cosmic string literature - multi¬ 
plied /i by Newton’s gravitational constant G in order to 
obtain a dimensionless number), a value which is below 
the current upper bound on G/x of G^ < 1.5 x 10“^ [1] 
(see also [S] for earlier work on limits on the string ten¬ 
sion). At this point, our string model had no more free 
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parameters. Interestingly, we found good agreement be¬ 
tween the predicted and observed globular cluster mass 
functions. 


However, in our previous study [T] we neglected the 
presence of string loop velocities (center of mass string 
loop velocities). Recent numerical simulations (eg. [5]) 
tell us that loops are typically born with translational 
velocities that are sizable fractions of the speed of light. 
The reason is that, since long string segments usually 
have relativistic speeds, then as string loops split off, 
the loops also gain significant velocities. When veloci¬ 
ties are taken into considerations, accretion will not be 
spherically symmetric. Additionally, accretion onto a 
moving loop may be less efficient compared to accre¬ 
tion onto a stationary loop. It should be noted, how¬ 
ever, that loop velocities also undergo red-shifting, and 
thus slow down as the loops age. 


In the following, we first review the hypothesis that 
globular clusters may be seeded by the cosmic string 
loops that arise from a string scaling solution in parti¬ 
cle physics models with a vacuum manifold which has 
the topology of a circle [T] . In section III we present our 
a first analysis of velocity effects on our model for the 
mass function of globular clusters. We compute a sup¬ 
pression factor which takes into account that loops with 
too large initial velocities will move a distance greater 
than the size of the spherical object which a station¬ 
ary loop would accrete, and we incorporate this factor 
into the predicted overall number density to determine 
a new mass function. This is then compared with the 
observed mass distribution of globular clusters in our 
Milky Way galaxy. Moving loops in fact do accrete 
matter, but do not give rise to a spherical distribution. 
In Section [Tv] we compare the total mass from spherical 
and non-spherical accretion, and determine an alterna¬ 
tive criterium for the maximal velocity of a loop that 
will seed a globular cluster. The resulting mass function 
turns out to be similar to the one derived using the first 
criterium. We then turn to a brief discussion of the cos¬ 
mic string rocket effect, another effect neglected in our 
previous work. We show that this effect has a negligi¬ 
ble impact on our globular cluster study. Finally, we 
present or conclusions in section |VI[ 
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II. GLOBULAR CLUSTERS FROM COSMIC 
STRING LOOPS 


In this report, we will consider a one-scale model for 
the distribution of strings which miH! implies that loops 
of initial radius Ri form at an initial time U given by 

ti = a~^l3Ri. ( 1 ) 

Here, a is a constant obtained from numerical solutions 
[ 6 ]. We will use the value a = 0.3. The average length of 
a string loop is given by I = /3R, where (3 is taken to be 
10 here; if the string loops were exactly circular /? would 
be 2tt. Making use of the fast decay approximation that 
assumes loops decay virtually instantaneously, we have 
R{t) = Ri{ti) until the decay time. String loops oscil¬ 
late and decay by the emission of gravitational radiation 
[9] whose strength is parametrized by a dimensionless 
constant 7 . Loops with radius smaller than 

Rc{t) = jGfit ( 2 ) 

live less than one Hubble expansion time before decay¬ 
ing. Hence, the number density of loops is constant for 
R < Rc{t). 

Linear cosmological perturbation theory tells us that 
accretion of matter around a cosmic string loop starts 
at teq, the time of equal matter and radiation. At this 
time, the mass function of the string loop with respect 
to radius is dominated by loops with a critical radius 
Rci , given by 

Rc^ = -fG^lteq ■ (3) 

We will use 7 = 10^. Though these cosmic string loops 
will have decayed by the present time, the objects they 
seed will continue to grow. The number density at a 
time t > teq for loops formed in the radiation dominated 
era of such objects inside a galaxy is given by (see e.g. 
m for reviews on the applications of cosmic strings in 
cosmology) 

n{R, t) = for R > R^^ 

n{R,t) = 

= const, for R < (4) 

The constant N depends on the square of the average 
number N of long string segments per Hubble volume 
since two long string segments are required to form a 
loop. We will take N = 10^. Incorporating the Zel- 
dovich approximation the local number density of 
string loops inside a galaxy will be enhanced by a fac¬ 
tor of F. This factor is estimated to be F = 64 due 
to accretion and virialization in each direction. In our 
calculations, we will use F = 10^. 

Assuming accretion continues to the present time, the 
mass which has accreted about these seed loops at the 


present time to is given by 

M(i?,,,to) = /37 (Gm)"^-/^'(^). (5) 

To obtain a feeling for the meaning of this expression, 
let us insert the values of to and G. We then have 

I ^ 10^^Mq , (6) 

where Mq stands for the solar mass. 

In our analysis, we take the peak mass (the mass 
where the observed globular cluster mass function for 
our Milky Way galaxy peaks) to fix our only free pa¬ 
rameter, the string tension. The mass function scales 
as for M > Me (which follows directly from 

the string loop distribution We predict a linear 

decay for M < Me- This comes about since the loop 
radius distribution is constant and loops with radius 
smaller than Re^ live only a fraction of a Hubble time 
step which scales linearly with R. 

In Fig. a comparison of the predicted mass function 
(the solid lines) with the observed distribution of glob¬ 
ular clusters in the Milky Way (histogram values) com¬ 
plied from [12] is made. To obtain the theoretical curve, 
we take the comoving number density n{R,teq)Yq 
loops (where n(i?, t) is given in multiply the result 
by the concentration factor F, allow each loop to grow 
in mass by a factor of Zeq (independent spherical accre¬ 
tion), and convert this into a mass distribution n{M, to), 
while taking into account the Jacobian of the transfor¬ 
mation from R to M. The result is then multiplied by 
the bin size 6M = fM, where / is a number, and by 
the volume V of the Milky Way galaxy. We obtain the 
following peak number density bin using <§ 

SN = ArF/aG2^-5/2^-3/2(g,^^-3/2^3/2^-3y^ ( 7 ) 

where z{t) is the cosmological redshift at time t. 

All of the calculations summarized in the this section 
assumed that cosmic strings loops are created and re¬ 
main at rest. In the following two sections, we will study 
the effect of velocities (translational center of mass mo¬ 
tion) through two different analyses. 

In our first analysis (discussed in the following sec¬ 
tion), we will compute the mass function of objects 
which accrete onto loops with velocities low enough such 
that the loop center moves a smaller physical distance 
than the physical radius of the shell which would be 
collapsing onto the loop if it had been stationary. If 
the loop moves further than this, we assume that no 
globular clusters forms. 

In our second analysis (to be discussed in the next 
to following section) we study the accretion of matter 
onto moving loops and keep only objects which are suf¬ 
ficiently spherical. We find that both conditions give 
similar resulting mass functions for globular clusters, 
the one from the second condition being slightly higher. 

In the following analyses, we are neglecting the rocket 
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FIG. 1. Dependence of the mass function of our model on the string tension G/r. The horizontal axis is mass on a logarithmic 
scale, the vertical axis gives the number density on a linear scale. The histogram show the data taken from [12]. The curves 
shown are for Gfj, = 2.92 x 10“^°, Gfx = 3.98 x 10“^°, Gfi = 5.43 x 10“^° and G/r = 7.41 x (in increasing order of 

mass at the peak position). The blue solid curve minimizes x^- The cosmic string parameters chosen are described in the 
text. The effects of string velocities are neglected. 


effect which is the effect of anisotropy in the loop gravi¬ 
tational radiation that causes loops to accelerate as they 
decay. We show in section |V] that this is indeed a good 
approximation. 


III. EFFECT OF COSMIC STRING 
VELOCITY: FIRST ANALYSIS 


of string loops for which 

AriR) < hiR,to), (10) 

and we assume that the accretion onto faster moving 
loops is not effective at producing compact globular 
clusters. Making use of ^ and we obtain an upper 
bound on the initial velocity: 


Accretion onto a cosmic string loop can be studied 
using the Zel’dovich approximation. As shown in e.g. 
in m, the physical distance h{q, to) from the center of a 
string loop to the mass shell which is “turning around” 
(i.e. becoming gravitationally bound) at the present 
time to is given by 

h(R.to) = (8) 

On the other hand, for an initial physical velocity Vi, the 
distance a loop of radius R has moved by the present 
time is 


Taking the distribution of initial velocities in each of 
three spatial directions to be a step function of width 
Vmax leads to the following probability that a string loop 
will satisfy the condition (11): 


V{v) = 


( 12 ) 


where Umax is the maximum velocity of cosmic strings 
determined through cosmic string evolution simula¬ 
tions. Taking the integral of the velocity distribution 
from zero to the upper bound on Vi, the rate of globular 
cluster formation becomes suppressed by a multiplica¬ 
tive factor S{R): 


Ar{R) = a{to) 


^ a{t) ) 


Vj 

liti) 


dt 


= 3a-i/2/3V2^i/4iy2^i/2^^ ^ 


S{R) = ia3/2^-l/2(g^)^l/4^y2^-3^^-l/2^ (13) 

The distribution has a cutoff when S{R) = 1 leading to 


In the present analysis we will only count the number 
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the critical radius i?c 2 • 


Rc2 = 


-6 


(14) 


Hence for values of i? < Rc^, there is no suppression 
from velocity effects and for R > Rc^ the suppression 
is given above by (13). Setting Rc^ = Rci we obtain a 


critical given by 


, (15) 

Thus, for z;max < 'f^max we have Rc^ > Ra and hence the 
mass function of predicted globular clusters from string 
loops will not change near the peak position compared 
to what was obtained in [1] neglecting the presence of 
string loop velocities. On the other hand, if Umax > '*^max 
then the mass function will be suppressed near the peak 
position. The relation of uJJj^x and Gp, is illustrated in 

Fig.i 

Taking into account the suppression factor when per¬ 
forming the calculations outlined in sectionjllj we obtain 
the following histogram of predicted number of globular 
clusters: 


For Rc^ > Rc^ : 

SN =NFfa^Gi3-5i2^-3i2 

X {Gfl)-^G^V2^-3y 

SN=i^)-^NFfa-^/3-^ 

X at Rc 2 

For Rc 2 < Rci ■ 

SN =^NFfa^r^j-^ 

X {Gfi)-^zlGtQ^v-l^V at i?ci 

SN ={^fNFfa^^Gi]-V2^-5/2 

X {G^l)-^GzlG^-3y-6^^Y at R^^ 


(16) 


(17) 


Notice that for i?ci < Rc 2 ^ the mass scales as M~^ 
for R > Rc 2 ■ In the radius R interval between i?ci and 
Rc 2 the mass function scales as M~^G as it does in the 
absence of velocity effects, and for masses smaller than 
Me, a linear decay is predicted by the same reasoning 
as in section [n] For R^^ > i?c 2 j the mass function scales 
as M“^ for R > as for R^^ < R < and 

decays linear for R < 77^2 ■ 

In Fig. we show that varying G/r shifts the peak 
position and amplitude of the mass for a fixed For 

Rc2 ^ Rci, initial velocity effects are negligible. How¬ 
ever, for Rc 2 < Rci but very close to the value of Rci 
there is a slight suppression in the region R^.^ < R < Rc^ 
from velocity effects. 

In Fig. 1^ we consider G/i = 5.43 x I0“^° which min¬ 
imizes in Fig- [ij we find from varying Umax that for 
^’max < 3.00 X 10“ , velocity has little effect on mass dis¬ 
tribution of globular clusters in the Milky Way galaxy. 
However, for Umax 3.00 x 10“^ we will not obtain a 
mass distribution. 
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FIG. 2. Analysis 1 - Relation between G/r and (15l. 
The horizontal axis is velocity on a linear scale and the ver¬ 
tical axis gives the G/r on a logarithmic scale. 


IV. EFFECT OF COSMIC STRING 
VELOCITY: SECOND ANALYSIS 


In the previous section we estimated the range of 
velocities for which spherical accretion onto a loop is 
a good approximation. On the other hand, accretion 
onto a moving loop can also be studied by means of the 
Zel’dovich approximation. This analysis is technically 
a bit more complicated than in the case of spherical ac¬ 
cretion, but the study has been carried out in |14] . The 
result is that half of the turnaround mass from a string 
with some initial velocity is within a region which can 
be approximated by a paraboloid of radius r — b^/^di 
and height h = 46^/^(ii), where (for loops born before 
the time of equal matter and radiation) 


b{t) 


I Gm 
15 Vlgteq 


a{t ), 


(18) 


where m is the mass of the loop. Note that since the 
accretion effectively starts at teq, it is the loop velocity 
Veq at that time which enters the formula. 

The value of the mass enclosed in this region is: 


= imait) for b{t) < 1. (19) 


Assuming that the accreted mass has uniform density, 
we find that density p is given by 

p = ^ = ^TT~^b~^d~^ma{t) for b{t) < 1. (20) 

Approximating the other half of the accreted mass to 
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Experimental • Theoretical 1 ' » ■ Theoretical 2 — • — Theoretical 3 —•—Theoretical 4 

FIG. 3. Analysis 1 - Dependence of the mass fnnction on G/i at Wmax = 3.00 x 10“^. The horizontal axis is mass on a 
logarithmic scale, the vertical axis gives the number density on a linear scale. The histogram shows data taken from m- 
The curves shown are for Gp = 2.92 x 10 Gfi = 3.98 x 10 G/i = 5.43 x 10 and Gfi = 7.41 x 10 (in increasing 
order of mass at the peak position). Notice that for the red and yellow solid curves, JJcj < Rci, for the blue solid curve 
Rc 2 = Rci and for the green solid curve Rc 2 > Rci- The blue solid curve minimizes for this particular Umax. 
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HMrt Experimental —•— Theoretical 1 —©—Theoretical 2 — • — Theoretical 3 • Theoretical 4 


FIG. 4. Analysis 1 - Dependence of the mass function at Gp = 5.43 x 10“^° on Umax. The curves are shown for Umax = 
1.00 X 10~^ (Theoretical 1), Umax = 6.46 x 10“^ (Theoretical 2), Umax = 3.44 x 10“^ (Theoretical 3) and Umax = 3.00 x 10“^ 
(Theoretical 4). Notice that for all Umax < 3.00 x 10“^ we would obtain the blue curve. The axes and data are the same as 
in the previous figure. 
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be spherical, the total mass is given by 


V. THE ROCKET EFFECT 


= ^ma{t). ( 21 ) 

Comparing this to the mass from spherical accretion: 

= ^ma{t) (22) 


we see non-spherical accretion results in a mass that is 
larger by a factor of two. 

Accretion is roughly spherical when the loop accre¬ 
tion sphericity parameter b(t) as defined in (181 is larger 
than one. In this analysis, we will consider a slightly 
lower bound by setting b{t) > 10“^. Using this condi¬ 
tion, we obtain an equation for v^q: 


V^q < (f )l/3/3l/3(G^)l/3^5/6^-l/3^1/3 ^ (23) 


The rocket effect is an effect which causes the string 
loop to accelerate as a consequence of anisotropic grav¬ 
itational radiation. Let us assume that the string loop 
has a center of mass velocity v^q at teq- Hubble expan¬ 
sion will lead to a decrease in the velocity whereas the 
anisotropic gravitational radiation will tend to cause an 
increase. Which effect wins depends on time. 

The equation of motion for the loop velocity taking 
into account Hubble expansion and anisotropic gravi¬ 
tational radiation is (see e.g. Appendix A of [TS] - we 
have generalized the result to times before teq) 


V{t) = Vi 


z{t) 

z{ti) 


5m 



tl/2 j 


, (26) 


Red-shifting the velocity to the time of loop formation 
in the radiation dominated era we obtain: 

U, < (|)l/3al/2/3-l/6(G^)l/3^1/12il/6^-l/6 ^ (24) 

Notice that the upper bound on the initial velocity in 
this analysis differs by only factor of 10^/^ from the 
upper bound found in the first analysis. From here, 
performing the same steps as in Analysis 1 would obtain 
results that are larger by a factor 10^/^. This can be 
seen clearly by determining the new 

. ( 25 ) 

The numerical results are presented in Fig. (§. 
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FIG. 5. Analysis 2 - Relation between and (251. 
The horizontal axis is velocity on a linear scale and the ver¬ 
tical axis gives the Gfi on a logarithmic scale. 


where Fp is the coefficient of anisotropic gravitational 
radiation, Mi is the mass at the loop formation time ti 
and m is the mass at a later time. We will take Fp ~ 
10. This value was estimated by numerical simulations 
using the Kibble-Turok loop solutions m- 

We now ask the question whether the rocket effect 
can have an important effect on gravitational accretion 
by the loops which dominate the globular cluster mass 
function. Thus, we are interested in loops with initial 
radius R = jCfiteq- If the second term on the right- 
hand side of (26) is smaller than the first term at the 


time t = teq, then we argue that the rocket effect will 
be negligible since even if at very late times the string 
loop starts to accelerate away from the initial position 
of the loop, the nonlinear object seeded by the loop is 
already in place one Hubble time after teq- Making use 
of (261 we then obtain the condition 


u,(Gm)-'/" > Tpa3/2p3/2^-3/2 ^ (27) 
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which, if satisfied, guarantees that the rocket effect is 
negligible. But for the values of Gn smaller than the 
current observational bound this condition is obviously 
satisfied. Hence, we conclude that the rocket effect does 
not effect our analysis. 


VI. CONCLUSIONS AND DISCUSSION 

In this paper we have considered the effects of string 
loop velocities on the model for the explanation of the 
origin of globular clusters proposed in ^ in which it 
is proposed that globular clusters are seeded by cosmic 
string loops. 

We see that velocities play a small, but noticeable role 
in the accretion of matter during the matter-dominated 
era. Cosmic string loops, born in the radiation era at ti 
with initial velocity Vi will travel a certain given by <§■ 
Demanding that this distance be smaller than the total 
size of the structure which accretes around a static loop, 
denoted in this paper by h{R,tf), leads to an upper 
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bound on the string velocities. In our Analysis 1 we 
assumed that loops with larger velocities do not give 
rise to globular clusters. We computed the corrections 
to the mass function obtained using this hypothesis, as 
a function of the maximal velocity of string loops. Our 
theory predicts a cutoff velocity. For velocities smaller 
than velocity effects are negligible. 

An improved analysis (Analysis 2) can be obtained 
by taking into account the non-spherical accretion of 
matter about a moving string loop, and demanding that 
the resulting nonlinear object is sufficiently spherical. 
The results of this analysis were very similar to those 
obtained using Analysis 1. 

Finally, we have shown that the rocket effect does not 


affect our analysis. 

Aside from globular clusters, it is noteworthy to men¬ 
tion that analyses discussed in this report may also be 
applied to ultra-compact-mini-halos, as they too can 
be manifestations of cosmic string loop accretion m 
Work on this topic is in progress. 
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